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SUMMARY 

This  document  is  a  report  of  the  measurements  of  the  light-scattering 
properties  of  blow-off  particulates  generated  In  simulated  thermonuclear  blasts. 
The  thermal  dust  particles  studied  were  generated  In  experiments  performed  by 
Science  Applications  Inc.  under  a  separate  DNA  contract*  using  a  solar  furnace 
In  Odelllo*  France.  The  samples  obtained  from  Science  Applications  Inc.  consist¬ 
ed  of  filters,  wall -wipings  and  loose  dust  taken  from -various  locations  In  the 
solar  furnace,  and  they  correspond  to  several  soil  samples  placed  In  the  furnace 
for  the  tests. 

The  data  reported  Include  particle  size,  complex  refractive  Index  and 
albedo  measurements  obtained  for  three  wavelengths  of  light  In  the  visible 
(blue-green)  range  of  the  spectrum  of  electromagnetic  radiation;  457.9,  488.0 
and  514.5  nm.  Earlier,  measurements  were  made  at  632.8  nm.  These 
light-scattering  data  were  obtained  by  suspending  single,  charged  particles  In 
an  electrodynamic  balance  In  the  path  of  a  polarized  laser  beam.  The  equipment, 
techniques  and  methods  of  data  analysis  are  outlined  herein. 

The  light-scattering  data  are  supplemented  by  scanning  electron  micrographs 
(SEM)  and  EOAX  (energy  dispersive  x-ray  spectrometry)  analyses,  which  provide 
Information  on  the  sizes,  shapes  and  chemical  composition  of  the  particles.  The 
EDAX  analyses  Indicate  that  most  of  the  dust  studied  was  mixed  silicates  of 
aluminum,  calcium.  Iron  and  potassium.  Particle  dimensions  ranged  from  about 
0.5  urn  to  100  Jim.  The  complex  refractive  Index  (m  =  m^  -  img)  was  found  to  have 
a  real  component,  m^.  In  the  range  1.481  -  1.723,  which  Is  consistent  with 
available  data  on  the  bulk  refractive  Index  of  silicates  and  clays.  The 
Imaginary  component,  m2,  was  small  (m2  <  0.06)  compared  with  results  available 
for  highly  absorbing  particulates  such  as  carbon  particles  (m2  =  0.6  -  1.0).  The 
particle  albedos  were  found  to  be  In  the  range  0.49  -  0.99. 
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PREFACE 


The  work  reported  here  represents  Phase  II  of  baseline  studies  on  the 
light-scattering  properties  of  blow-off  particulates.  Phase  I  Involved  the 
development  and  testing  of  equipment,  techniques  and  computer  software,  and  the 
results  of  Phase  I  were  summarized  In  Report  number  DNA  6274F,  dated  April  1983. 
That  Initial  work  was  supported  under  Contract  Number  DNA  001-80-C-0213. 

Phase  II  was  continued  at  the  Institute  of  Paper  Chemistry  (IPC)  In 
Appleton,  Wisconsin,  but  In  August  1983  the  Principal  Investigator,  Dr.  E.  James 
Davis,  joined  the  faculty  of  chemical  engineering  at  the  University  of 
Washington.  His  light-scattering  laboratory  was  re-established  at  the  Universi¬ 
ty,  and  the  research  on  Phase  II  was  continued  there  under  a  subcontract  with 
IPC. 

This  research  provides  Information  on  the  scattering  and  absorption  prop¬ 
erties  In  the  visible  region  of  the  spectrum  of  electromagnetic  radiation  needed 
to  Interpret  and  assess  field  measurements  of  dust  particles  generated  In  a 
simulated  thermonuclear  explosion. 
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SECTION  1 


OBJECTIVES 

The  research  objectives  were  to  measure  the  light-scattering  properties  of 
dust  particles  generated  thermally  in  solar  furnace  experiments  performed  by 
Science  Applications  Inc.,  and  to  determine  the  effects  of  particle  size  and 
shape  and  soil  sample  on  the  optical  properties.  The  properties  determined  are: 

1.  Phase  function  (intensity  vs  scattering  angle)  for  dust  particles  in 

o  o 

the  wavelength  range  3500  A  -  7000  A. 

2.  Complex  refractive  index. 

3.  Particle  albedo  (or  absorption  and  extinction)  for  various  wave¬ 
lengths. 

1.1  BACKGROUND 

1.1.1  Light  Scattering  Theory 

The  optical  properties  of  the  blow-off  particulates  were  determined  from 
light  scattering  data  using  Mie  theory  to  interpret  the  data.  From  Mie  theory 
the  intensity  of  the  scattered  light  is  established  as  a  function  of  sphere 
radius,  a,  wavelength,  x,  complex  refractive  index,  m  =  m^  -  img,  and  scattering 
angle,  e.  The  dimensionless  size  of  the  sphere,  a,  is  defined  by  a  =  2ira/x.  If  a 
and  m  are  known,  the  scattered  intensity,  1^,  as  a  function  of  angle  e  for  light 
polarized  with  the  electric  vector  normal  to  the  plane  of  scattering  (the  TM 
mode)  are  calculated  from  the  equation. 


where 


h  ’  TiTHTIT  ®>  *  ‘‘n  'n  • 


ir^(cos  e)  =  P^^^(cos  0)/sin  e 


T„(cos  e)  =  pjj^^cos  e) 


where 


D^(b)  -  m  D^(o)  jj^(o) 

‘n  ~  0  (b)  -  m  G  (o)  h  (o) 
n  n  n 


m  D^(b)  -  D^(o)  j^(o) 

m  D  (b)  “  G  (o)  h  (a) 
n  n  n 


h^(o)  =  j^(o)  +  i  N^(a) 


s„(2)  -  (zh„)Vzh„ 


In  the  above  equations  h^  Is  the  Hankel  function,  Is  the  spherical  Bessel 
function,  Is  the  Neumann  function,  P^^^  Is  the  Legendre  function,  and  6  *  ma. 

Figures  1  and  2  show  results  of  Mie  theory  calculations  as  three  dimension¬ 
al  plots  of  the  Intensity  ratio  Ij/Iq  {where  Is  the  Incident  beam  Intensity) 


as  functions  of  the  other  parameters.  Figure  1  demonstrates  the  effects  of 
particle  size  and  angle  on  for  a  real  refractive  index  m^  =  1.696  (m2  =  0  for 
Fig.  1),  and  Fig.  2  shows  the  effects  of  the  imaginary  component,  m2,  and  angle 
for  mj  =  1.550  and  a  =  10.  Figure  1  indicates  that  as  o  increases  the  number  and 
intensity  of  the  peaks  in  the  angle  range  shown  increase.  Furthermore,  the 
angular  positions  of  the  peaks  and  troughs  depend  on  o  (and  m).  This  is  what 
makes  it  possible  to  determine  the  size  and  refractive  index  from 
light-scattering  data  on  versus  e. 


Figure  1.  The  dependence  of  scattered  intensity  1^,  on  optical  size  a 
and  angle  0  for  a  refractive  index  m  ®  1.696  -  iO. 


Figure  2  illustrates  an  important  consideration  in  the  evaluation  of  m,.  As 


mg  increases  from  zero  (no  energy  absorption),  the  intensity  profile  becomes 
nearly  independent  of  mg,  and  it  becomes  difficult  to  determine  mg  (and,  hence, 


Intensity  for  m^  =  1.55  and  a  =  10. 


The  complicated  topology  shown  in  Figures  1  and  2  must  be  taken  into 


account  in  the  analysis  and  interpretation  of  the  data,  for  since  a  and  m  are 
Initially  unknown,  the  Mie  theory  equations  must  be  used  In  an  Iterative  scheme 
to  determine  them  from  the  measured  scattering  function  I^(e).  The  data  analysis 
procedure  Is  outlined  below. 

1.1.2  Principles  of  the  Electrodynamic  Balance 

The  scattering  Intensity  profile,  (e),  has  been  measured  by  suspending  a 
single,  charged  dust  particle  In  an  electrodynamic  balance/light-scattering 
photometer  developed  by  Davis  and  Ray^.  A  cross  section  of  the  electrodynamic 
balance  Is  shown  In  Figure  3,  and  an  exploded  view  of  the  device  Is  shown  In 
Figure  4.  The  endcap  electrodes  have  geometries  described  by  the  equation 


where  Zz^  Is  the  minimum  distance  between  the  endcap  and  the  ring  electrode 
geometry  satisfies  the  equation 


The  ring  electrode  was  connected  to  an  ac  power  supply  (usually  a  60  hz 
transformer)  to  focus  the  charged  particulate  on  the  centerline  (vertical  axis) 
of  the  balance.  With  no  dc  field  applied  the  particle  oscillates  about  a  point 
below  the  midplane  of  the  balance.  When  a  dc  field  Is  applied  to  offset  the 
gravitational  force  on  the  particle  It  remains  stationary  at  the  midplane  (z=0) 
between  the  endcaps. 


Davis,  E.J.  and  Ray,  A.K. ,  "Single  Aerosol  Particle  Size  and  Mass 
Measurements  Using  an  Electrodynamic  Balance,"  0.  Colloid  and  Interface 
Scl.,  75(2),  566-576,  1980. 
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Figure  3.  The  electrodynamic  balance/light  scattering  photometer 
and  peripheral  equipment. 
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Figure  4.  An  exploded  view  of  the  electrodynamic  balance 


If,  for  a  negatively  charged  particle,  a  dc  voltage,  is  applied  to  the 
upper  electrode  and  is  applied  to  the  lower  endcap,  the  electric  field  near 
the  null  point  (midplane)  of  the  balance  is  given  by 


where  C  is  a  geometrical  constant  which  takes  into  account  the  shape  of  the 
electrodes  and  the  edge  effects. 

For  a  particle  of  mas5,  m^,  and  charge,  q,  the  particle  is  balanced  when 
the  following  force  balance  is  satisfied 


“lEdc.z  -  ”p9  -  0. 


where  g  is  the  gravitational  acceleration  constant.  The  dc  voltage  required  to 
balance  the  particle  is  given  by 


u  -  P  0 
dc  qC 


As  indicated  in  Figure  4,  a  laser  beam  was  aligned  to  pass  through  the 
center  of  the  chamber,  and  a  window  and  slit  in  the  ring  electrode  permitted 
detection  of  the  scattered  light  by  means  of  a  transversing  photomultiplier 
tube.  The  photodetector  was  driven  by  a  stepping  motor,  and  the  signals  from  the 
photodetector  (intensity)  and  the  stepping  motor  (angle)  were  digitized  and 
recorded  using  a  Digital  Equipment  Corporation  MINC-23  computer.  Additional 
details  of  the  experimental  system  are  provided  in  the  next  section. 


SECTION  2 


THE  EXPERIMENTAL  APPARATUS  AND  DATA  ANALYSIS 

2.1  THE  EXPERIMENTAL  SYSTEM 

The  components  of  the  experimental  system  are  indicated  in  Figure  5.  The  dc 
power  supplies  were  matched  and  connected  to  a  dual  potentiometer  to  permit 
precise  balancing  of  the  charged  particulate,  observed  through  a  low-power 
microscope  attached  to  the  balance  frame.  Dust  particles  were  washed  from  the 
filters  or  wall-wiping  material  to  form  a  dilute  suspension  in  water,  and  the 
suspension  WbS  drawn  into  a  microliter  syringe.  The  needle  of  the  syringe  was 
connected  to  a  high  voltage  dc  supply  to  facilitate  charging. 


Injection  of  a  dust  particle  into  the  balance  was  accomplished  by  forming  a 


small  droplet  at  the  top  of  the  needle,  switching  on  the  high  voltage,  and 
adjusting  the  ac  voltage  to  "catch"  the  droplet.  The  water  surrounding  the 
particulate  quickly  evaporated,  leaving  the  charge  on  the  particle.  The  dc 
voltage  was  then  applied  to  the  endcaps  to  stabilize  the  particle  at  the  center 
of  the  balance.  Usually  a  dust  particle  was  suspended  for  more  than  one  hour 
prior  to  obtaining  a  light-scattering  data  to  insure  that  the  water  was  removed. 

After  no  mass  change  could  be  detected,  the  laser  intensity  and  voltage 
input  to  the  photomultiplier  tube  were  adjusted  to  give  a  favorable  sig¬ 
nal  -to-noise  ratio,  indicated  by  recording  the  scattered  intensity  as  a  function 
of  angle  on  an  X-Y  recorder.  Next  the  computerized  data  acquisition  was  carried 
out. 

The  laser  light  sources  provided  vertically  polarized  (TM  mode)  light  at 
any  of  several  discrete  wavelengths.  Two  sources  were  available:  (i)  a  10  mW 
Spectra  Physics  helium-neon  laser  (632.8  nm)  and  (ii)  a  Coherent  Innova  90-5  5W 
argon-ion  laser.  A  tunable  dye  laser  was  available,  but  was  not  used  for  the 
results  reported  here  because  the  argon-ion  laser  had  a  sufficient  range  of 
wavelengths  in  the  visible.  The  wavelengths  used  were  457.9,  488.0  and  514.5  nm. 

For  data  acquisition  the  signals  from  the  photomultiplier  tube  and  the 
stepping  motor  were  conditioned  to  be  in  the  voltage  range  0-5  V  by  means  of 
operational  amplifiers,  and  these  signals  were  digitized  and  recorded  using  a 
DEC  MINC-23  computer  system.  The  computer  data  acquisition  was  activated  by 
triggering  the  computer  clock  with  a  signal  from  the  stepping  motor,  and  the 
sampling  rate  was  selected  to  obtain  about  500  data  points  during  a  forward  scan 
(from  about  6®  to  174®)  and  the  same  number  during  a  reverse  scan.  Data-taking 
stopped  after  the  completion  of  one  forward-reverse  scan  cycle. 

Figure  6  shows  typical  raw  data  for  a  forward  scan,  plotted  on  the 
HP  7245  B  Printer/Plotter  as  intensity  (PMT  output  voltage)  versus  angle 


(stepping  motor  voltage).  Two  sharp  troughs  should  be  noted  at  about  1.35  volts 
and  3.78  volts.  These  minima  in  the  scattered  intensity  were  caused  by  the 
blockage  of  the  light  by  two  pegs  installed  to  provide  calibration  of  the 
voltage  signal  from  the  stepping  motor.  The  raw  data  (PUT  output  voltage  versus 
stepping  motor  voltage)  were  examined,  using  appropriate  computer  software,  to 
determine  the  voltages  at  which  the  sharp  troughs  occurred.  These  stepping  motor 
voltages,  then,  corresponded  to  the  angles  16.2“  and  163.8“. 


RAW  LIGHT  SCATTERING  DATA  FOR  RUN  SB  RB5I-2  WALL  COLLECTION 
wavelength  •  457.  B  nm 
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Figure  6.  Raw  light  scattering  data  from  a  dust  particle,  A  =  457.9  nm 


After  a  scan  cycle  was  completed  the  wavelength  of  the  laser  was  changed, 
and  the  experiment  repeated.  At  the  end  of  a  sequence  of  experiments  on  a 
particular  dust  particle  the  particle  was  expelled  from  the  balance,  and  the 


background  signal  was  measured  and  recorded  at  each  wavelength.  Figure  6  shows 
both  the  data  obtained  with  a  dust  particle  suspended  in  the  beam  and  the 
background.  The  data  used  in  further  analysis  was  the  difference  between  the  two 
signals. 

2.2  INTERPRETATION  OF  THE  DATA 


To  determine  the  particle  size  and  refractive  index  from  data  of  the  type 
shown  in  Figure  5  Mie  theory  was  applied.  Forty  data  points  uniformly  spaced  on 
the  angle  axis  were  selected,  usually  in  the  range  30®  -  130°,  and  the  inten¬ 
sities,  I^  predicted  from  Mie  theory  for  these  angles  were  calculated.  The 
voltages  from  the  photodetector  circuit  were  converted  to  normalized  intensities 
using  the  relation 

le(0i)  *  o[V{e^.)  -  Vj,(0.)],  (14) 

where  o  is  a  scale  factor,  V{0^. )  is  the  signal  at  the  i^^  angular  position, 
Vb(0i)  is  the  background  signal  (no  particle  present)  at  that  position  and 
lg(0^)  is  the  experimental  intensity  at  e^.  The  scale  factor,  a,  can  be  de¬ 
termined  in  a  number  of  different  ways.  One  method  is  to  use  the  maximum  (high¬ 
est  peak  in  the  light  scattering  profile)  and  minimum  (lowest  trough  in  the 
profile)  detector  signals  to  normalize  lg(0^)»  that  is,  define  a  by 


a  - 


ly 


1 


max 


mm' 


(15) 


This  method  places  too  much  emphasis  on  only  two  points  of  the  light 
scattering  profile,  for  if  there  is  excessive  noise  in  the  signal,  significant 


error  is  Introduced  in  the  normalization.  A  better  approach  is  to  consider  a  a 
parameter  in  the  optimization  procedure  used  to  compare  the  experimental  inten¬ 
sities  with  those  computed  from  Mie  theory. 


To  interpret  the  experimental  data  we  determined  values  of  o,  m 
(m  =  m^  -  img)  and  o  that  minimize  the  function 

N 

F(a,m,o)  =  I  w. 

i=l  ^ 

where  f^.  is  the  i^*^  residual  defined  by 


(16) 


fi  =  Ie(e,)  - 


(17) 


and  w^.  is  the  weighting  factor.  For  these  studies  N  was  usually  40,  and  the  most 

2 

suitable  weighting  factor  was  found  to  be  w.  *  sin  e^. 

To  minimize  the  objective  function,  F(a,m,a)  we  applied  the 
Levenberg-Marquardt  algorithm  (Marquardt  ),  which  performs  an  optimum  inter¬ 
polation  between  the  Taylor  series  and  gradient  (steepest  descent)  methods.  The 
Gauss-Newton  method  of  optimization  depends  upon  reduction  of  the  residual,  f^, 
to  linear  form  by  first  order  Taylor  approximations  taken  about  an  initial  or 
trial  solution  for  the  parameters.  This  method  frequently  fails  because  of 
divergence  of  successive  iterations,  so  various  modifications  of  the  method  of 
steepest  descent  have  been  introduced  to  avoid  such  divergence,  but  the  method 
of  steepest  descent  can  give  agonizingly  slow  convergence  after  the  first  few 


Marquardt,  D.W. ,  "An  Algorithm  for  Least-Squares  Estimation  of 
Nonlinear  Parameters,"  J.  of  the  Society  of  Industrial  Applied  Mathematics, 
11(2),  431-441,  1963. 
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RESULTS 


3.1  SAMPLE  INDENTIFICATION 


To  identify  the  dust  samples  in  this  report  we  shall  use  the  run  numbers 
used  by  Science  Applications,  Inc.  [Trip  Reports,  March  1980  and  September  1980, 
DNA  Soil  Test  Program  at  the  CNRS  1000  KW  Solar  Furnace,  Odeillo,  France, 
November  24,  1980].  Table  1  lists  the  run  numbers,  soil  type  codes  and  descrip¬ 
tion  for  each  of  the  samples  studied  herein. 

The  samples  described  above  as  wall  wipings  were  particulates  washed  from 
towelling  used  to  remove  dust  from  the  furnace  walls.  Very  little  dust  was  found 
in  filter  samples,  and  it  was  difficult  to  wash  dust  out  of  the  filter  matrix. 


3.2  ELECTRON  MICROSCOPY 

A  scanning  electron  microscope  equipped  for  elemental  analysis  by  energy 
dispersion  x-ray  spectrometry  (EDAX)  was  used  to  examine  the  sizes,  shapes  and 
elemental  analyses  of  dust  particles  in  the  samples.  Figures  7  and  8  are  scan¬ 
ning  electron  micrographs  (SEM)  of  loose  thermal  dust  collected  from  the  solar 
furnace  and  wall  wipings,  respectively,  for  Run  57,  and  they  illustrate  the 
range  of  sizes  (note  the  100  ym  bars  at  the  bottom  of  each  picture)  and  the 
shapes  encountered.  Approximately  half  of  the  particulates  shown  in  these 


I 


figures  are  spherical  or  near-spherical,  indicating  that  the  temperatures 
reached  were  sufficient  to  melt  many  of  the  chemical  species  present. 


Table  1.  Identification  of  samples  studied 


Run  No. 

Soil  Type  Code 

Description 

57 

Silt  P-14 

Wall  wipings  and 

Bottom  filter 

60 

Clay  P-13 

Wall  wipings 

69 

RB5I-2 

Wall  wipings  and 

Bottom  filter 

101 

Wall  wipings 

104 

RC5Y-2 

Wall  wipings 

107 

RC4U-2 

Wall  wipings 

108 

RD5I-2 

Wall  wipings 

136 

Soil  RBU2#15 

Wall  wipings 

138 

Gray  tile 

Wall  wipings 

147 

Local  vegetation 

Wall  wipings  and 

Bottom  filter 

1180 

5F 

Wall  wipings 
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Figures  9  and  10,  on  the  otherhand,  show  many  particles  of  irregular  shape 
and  relatively  large  size  (dimensions  >  100  ym)  for  loose  material  from  Runs  69 
and  1180,  respectively.  Such  large  particles  would  be  expected  to  have  fairly 
rapid  settling  rates  and  were  not  found  in  the  filters.  It  is  clear  from  Fig¬ 
ure  7-10  that  the  range  of  particle  sizes  is  very  great,  and  Figure  11-13  for 
Run  136  show  that  wide  ranges  of  sizes  and  shapes  are  encountered  in  a  single 
run.  The  three  figures  were  obtained  by  increasing  the  magnification  from  X200 
to  X2400.  At  a  magnification  of  X200  many  agglomerated  particulates  can  be  seen 
(Figure  11),  and  at  the  higher  magnifications  a  very  large  number  of  irregularly 
shaped  particles  are  seen. 

A  set  of  SEM  and  the  corresponding  EDAX  analyses  are  shown  in  Figure  14 
through  26  for  the  runs  listed  in  Table  1.  Several  of  the  EDAX  analyses  show 
large  peaks  corresponding  to  gold  (Au)  and  palladium  (Pd),  for  these  elements 
were  used  to  shade  the  particulates  for  some  of  the  SEM.  The  SEM  selected  are 
for  various  magnifications;  some  show  relatively  large  regions  of  a  sample,  and 
others  focus  on  specific  particulates.  Table  2  summarizes  the  main  elements 
detected  in  the  samples. 
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Table  2.  Elemental  analysis  of  various  samples 
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Run  No. 


Main  Elements 


A1,Ca,K.Na,S1 
A1 ,Ca,Fe,K,S1 
A1,Ca.K.Na,S1 
A1  .SI 
Al.SI 
Al.Ca.SI 


Al.Ca.Fe.K 
Ag.AI .Ca.SI 
Al.Cl.K.Mg.SI 
A1  ,Ca,S1 


Trace  Elements 


Ag.Fe.Cu 


Cu  .Fe 


Most  of  the  samples  consisted  of  mixed  silicates  of  aluminum,  calcium  and 
potassium,  and  samples  from  Runs  57  and  69  contained  sodium  compounds.  Iron  was 
found  In  significant  amounts  In  samples  from  Runs  60  and  136,  and  Iron  and 
copper  were  present  In  trace  amounts  In  most  samples.  The  Iron  and  copper  might 
have  been  in  the  form  of  oxides,  but  quantitative  analysis  and  Information  on 
the  chemical  species  encountered  were  not  possible  with  the  techniques  used. 

Samples  from  Run  147  were  found  to  be  different  from  the  others  in  that 
magnesium  and  chlorine  were  detected  In  these  samples  which  came  from  local 
vegetation.  Run  136  has  some  unique  features,  for  Figure  11-13  show  agglomerated 
particles  and  particles  that  might  be  of  biological  origin  (e.g.,  spaces). 
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Because  spherical  particulates  are  not  seen  (see  Figure  13  and  23a)  it  appears 
that  the  temperatures  encountered  in  Run  136  were  too  low  to  melt  silicates. 


3.3  LIGHT  SCATTERING  PHASE  FUNCTIONS 


Figure  27-32  are  representative  plots  of  scattered  intensity  versus  angle 
for  various  samples  at  various  wavelengths,  and  a  complete  set  of  raw  data  are 
given  in  Appendix  A.  The  figures  shown  in  this  section  correspond  to  data  for 
which  the  optimization  program  converged  to  yield  values  of  size  a  =  irDp/A  and 
refractive  index  m  =  m^  -  im2.  For  many  experiments  no  satisfactory  match 
between  the  measured  phase  functions  and  Mie  theory  was  achieved,  so  for  experi¬ 
ments  only  raw  data  are  given  in  the  appendix. 

For  many  runs  only  the  angle  range  from  about  20°  to  100°  was  examined,  but 
with  the  redesign  of  the  viewing  window  more  recent  data  cover  the  range  from 
about  20°  to  160°.  The  pegs  used  to  calibrate  the  stepping  motor  voltage  pre¬ 
vented  the  collection  of  data  beyond  this  range. 

Figures  27-29  show  a  set  of  data  for  a  thermal  dust  particle  from  Run  60 
obtained  at  three  different  wavelengths.  The  agreement  between  experimental  and 
calculated  phase  functions  is  fairly  good  here,  considering  that  the  data  are 
"noisy",  but  there  is  a  rather  larger  scatter  in  the  values  of  particle  diame¬ 
ter,  Dp,  calculated  from  the  values  of  a  obtained  from  the  data  analysis.  The 
calculated  particle  diameters  are  0.730  pm,  0.763  pm  and  0.881  pm  with  a  mean  of 
0.791  pm.  The  scatter  about  the  mean  is  about  ±10%.  There  is  also  scatter  in  the 
values  of  m^,  and  m2  from  the  data,  but  the  refractive  index  can  be  expected  to 
vary  with  wavelength,  and  more  will  be  said  on  that  point  below. 

A  second  set  of  three  runs,  which  cover  a  different  range  of  angles,  is 
presented  as  Figures  30-32  for  Run  1180.  Note  that  there  is  backscatter  shown  in 
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these  figures.  The  particle  diameters  calculated  for  this  particle  from  the 
three  experiments  are  1.401  um,  1.429  um  and  1.523  urn  with  a  mean  of  1.451  urn. 
the  scatter  is  only  about  ±5%  for  this  set.  There  is  rather  little  variation  in 
refractive  index  for  these  three  runs,  as  shown  in  Table  3,  which  summarizes  the 
results  for  the  runs  for  which  convergence  of  the  matching  procedure  occurred. 
Also  given  in  the  table  are  the  calculated  albedos  for  each  run  and  the  average 
values  of  a,  m^,  m2  for  each  particle  measured  at  more  than  one  wavelength. 
Average  values  of  the  parameters  permit  us  to  compute  more  comprehensive  phase 
function  plots,  and  Figures  33-35  show  composite  phase  functions  for  particles 
from  several  runs.  These  figures  are  based  on  the  mean  values  of  a,  and  m2 
reported  in  Table  3  for  each  set  of  data  on  a  particular  dust  particle. 

Additional  comparisons  between  measured  phase  functions  and  Mie  theory 
results  are  given  in  the  appendix,  and  many  of  the  comparisons  show  poor  agree¬ 
ment  with  the  detailed  structure  of  the  intensity  profile  but  broad  agreement 
with  the  positions  of  peaks,  troughs  and  general  structure.  This  might  be 
expected  of  particulates  that  are  not  perfect  spheres  but  have  blemishes  due  to 
outgassing  or  other  phenomena  that  occurred  during  solidification  after  the 
particle  melted. 


Figure  30.  Experimental  phase  function  data  (*)  for  a  particle  from  Run 
1180  for  A  =  457.9  nm  compared  with  Mie  theory  ( — )  for  a  = 
and  m  =  1.719  -  i  0.029. 
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SECTION  4 


DISCUSSION  OF  RESULTS 

Table  3  shows  that  the  values  of  determined  by  comparing  phase  function 
data  with  Mie  theory  are  relatively  small  (compared  with  strongly  absorbing 
particulates  such  as  carbon  or  soot  particles  for  which  m2  >  0.6)  and  scattered. 
Since  mg  strongly  affects  the  computed  albedo*  the  latter  Is  also  scattered.  For 
the  small  values  of  mg  encountered  here  It  Is  difficult  to  determine  mg  from 
noisy  phase  function  data  for  the  reasons  suggested  by  Figure  36,  which  shows 
calculated  phase  functions  for  a  sphere  with  a  -  10,  m^  =  1.600  and  various 
values  of  mg.  The  primary  effect  of  Increasing  mg  Is  to  decrease  the  Intensity 
of  the  scattered  light  without  appreciably  shifting  the  angular  positions  of 
peaks  and  troughs.  Thus,  It  Is  difficult  to  extract  accurate  values  of  mg  by 
Inverting  the  phase  function  Information.  The  values  of  mg  reported  here  should 
be  construed  to  be  only  order  of  magnitude  estimates. 

The  range  of  dust  particle  sizes  studied  with  the  electrodynamic  balance  Is 
considerably  smaller  than  the  range  of  particles  shown  in  the  scanning  electron 
micrographs.  There  are  two  likely  reasons  for  this  discrepancy  (1)  the  number 
density  of  micron-size  particulates  Is  very  much  greater  than  that  of  larger 
particulates  (say,  Dp  >  20  wm),  and  (11)  larger  particulates  Injected  Into  the 
electrodynamic  balance  have  much  larger  settling  velocities  than  small  ones,  so 
the  latter  are  preferentially  captured.  Efforts  were  made  to  randomize  the 
sampling,  but  few  large  particulates  were  encountered  In  the  light-scattering 
studies. 

More  than  one-half  of  the  particulates  studied  In  the  electrodynamic 
balance  could  not  be  compared  with  Mie  theory,  for  they  were  not  spherical.  One 
Indication  of  asymmetry  of  a  particulate  was  that  when  observed  through  the 


Table  3.  Results  for  dust  particles  for  three  wavelengths  (Continued) 


RUN 

X,nm 

a 

mi 

m2 

A1 bedo 

a  ,um 

Run  104  RC5Y-2  Wall  wipings 

457.9 

12.23 

1.601 

0.024 

0.5589 

0.8915 

=  888.0 
ac 

488.0 

12.18 

1.593 

0.007 

0.7619 

0.9463 

=  0.906 

514.5 

12.38 

1.573 

0.008 

0.7243 

1.014 

AVERAGE 

12.26 

1.589 

0.013 

0.6817 

0.9506 

Run  107  P-14  Wall  wiper 

sample  RC4U-3 

457.9 

13.65* 

1.594* 

0.006* 

0.8028* 

0.9948* 

Vac  =  768.0 

488.0 

12.40 

1.574 

0.007 

0.7417 

0.9631 

Vdc=  0.792 

514.5 

12.53 

1.559  . 

0.005 

0.7816 

1.026 

AVERAGE 

12.86 

1.576 

0.006 

0.7754 

0.9946 

Run  108  Wall  wipings 

457.9 

12.75 

1.529 

0.005 

0.8013 

0.9291 

Vac  =744.0 

488.0 

12.88 

1.510 

0.005 

0.796 

1.000 

Vdc  =  0.675 

514.5 

13.08 

1.481 

0.007 

0.78 

1.0711 

AVERAGE 

12.90 

1.507 

0.006 

0.7924 

1.000 

Run  136  RBU2  #15  Wall 

wipings 

457.9 

- 

- 

- 

- 

- 

V^^  =  1032.0 

488.0 

12.40 

1.582 

0.002 

0.8861 

0.9632 

Vjc  =  0.695 

514.5 

- 

- 

- 

- 

- 

AVERAGE 

12.40 

1.582 

0.002 

0.8861 

0.9632 

*  extrapolated 
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Table  3.  Results  for  dust  particles  for  three  wavelengths  (Concluded) 


RUN 

X,nm 

a 

mi 

Al bedo 

a  ,ym 

Run  138  Black  Tile  Wall 

wipings 

457.9 

14.14 

1.723 

0.007 

0.7091 

1.0303 

Vac  = 

488.0 

14.59 

1.674 

0.010 

0.6716 

1.1332 

^dc'  1-27 

514.5 

15.06 

1.613 

0.000 

0.9892 

1.2333 

AVERAGE 

14.60 

1.670 

0.006 

0.79 

1.1323 

Run  1180  Wall  5-F 

457.9 

9.61 

1.719 

0.029 

0.5531 

0.7004 

Vac  =  1824.0 

488.0 

9.20 

1.716 

0.001 

0.9369 

0.7145 

Vdc=  82.0 

514.5 

9.30 

1.709 

0.026 

0.5768 

0.7615 

AVERAGE 

9.37 

1.715 

0.019 

0.689 

0.7255 

microscope  eyepiece  on  the  picobalance  the  particle  ''twinkled''  due  to  irregular 
scattering  as  it  rotated.  Light-scattering  data  taken  with  such  particles  showed 
no  reproducibility  during  forward  and  backward  scans  of  the  traversing  photo¬ 
multiplier  tube,  and  those  data  are  not  included  in  this  report. 

The  real  parts  of  the  refractive  indices  determined  in  this  research  are 
generally  consistent  with  refractive  indices  reported  for  minerals  (silicates, 
clay,  mica,  etc.).  Table  4  lists  typical  values  of  m^  for  feldspars,  clay,  mica 
and  other  silicates  containing  Al,  Ca,  Mg  and  Na.  The  range  of  values  listed  in 
Table  4  is  about  the  same  as  that  determined  here  for  the  various  soil  samples. 

There  is  no  clearcut  dependence  of  refractive  index  on  wavelength  for  the 


wavelengths  studied  here  and  in  Phase  I  of  the  project  (457.9  nm  -  632.8  nm), 
but  is  is  to  be  expected  that  in  the  ultraviolet  region  of  the  spectrum  is 
larger,  for  the  silicates  can  be  expected  to  absorb  in  the  UV  range. 


There  is  also  no  well-defined  variation  in  optical  properties  from  one  soil 
sample  to  another,  for  variations  between  individual  particles  in  a  single 
sample  are  large.  That  is  clear  from  an  examination  of  the  SEM  data  in  Sec¬ 
tion  3,  for  the  particulates  vary  greatly  in  size,  structure  and  melting  points. 

Table  4.  Refractive  indices  of  minerals 


NAME 

FORMULA 

REFRACTIVE  INDEX,  m 

Andes ine 

(Ca0,Na20)Ai203*4Si02) 

1.549 

Anorthite 

Ca0*Ai203*2Si02 

1.5755 

Anorthoclase 

(Na,K)20-Ai203-6Si02 

1.523 

Laboradorite 

NaAiSi30g  +  CaAi2Si20g 

1.559 

Oligoclase 

NaAiSigOg  +  CaAigSiOg 

1.539 

Orthoclase 

K20.Ai203»6Si02 

1.518 

Wollastonite 

CaSiOg 

1.616 

Kaolinite  (china  clay) 

Ai203*2Si02*2H20 

1.561 

Muscovite  (white  mica) 

K20-3Ai203-6Si02-2H20 

1.561 

Anal  cite 

NaAlSi20g-H20 

1.479-1.493 

Andalusite 

Ai20Si0^ 

1.629-1.640 

Clinozoisite 

Ca2Al3Si30j2(0H) 

1.670-1.715 

Enstatite 

MgSiOg 

1.650-1.662 

Forsterite 

Mg2Si04 

1.635 

Gehlenite 

Ca2Al2Si0j 

1.669 

Kyanite 

Al20Si0^ 

1.712-1.718 
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APPENDIX  I 

EXPERIMENTAL  PHASE  FUNCTION  DIAGRAMS  FOR  THE  DUST  PARTICLES 


The  following  phase  function  diagrams  are  presented  in  the  same  order  of 
increasing  run  number  used  for  Figures  14-26.  In  some  cases  the  experimental 
data  are  compared  with  Mie  theory,  but  in  others  the  inversion  of  the  phase 
function  data  by  means  of  the  iterative  optimization  procedure  discussed  in 
Section  2  did  not  converge,  so  the  parameters  a,  m^  and  m2  could  not  be  de¬ 
termined. 
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Figure  Al.  Experimental  phase  function  data  (*)  for  a  particle  from 
57  for  X  =  4  57.9  nm  compared  with  Mie  theory  ( — )  for  a  ■■ 
and  m  =  1.55  -  i0.003. 


Experimental  phase  function  data  (*)  for 
57  for  X  =  457.9  nm. 


ANGLE  IN  DECREES 

Experimental  phase  function  data  (*)  for  a  particle  from  Run 


Experimental  phase  function  data  {*)  for  a  particle  from  Run 
60  for  X  =  514.5  nm  compared  with  Mie  theory  ( — )  for  oi  =  4 . 
and  m  =  1.633  -  i0.026. 


Experimental  phase  function  data  (*)  for  a  particle  from  Run 


601  -  i  0.024 


Experimental  phase  function  data  (*)  for  a  particle  from  Run 


Figure  ^^31.  Experimental  phase  function  data  (*)  for  a  particle  from  Pun 
138  for  ^  =  457.9  nm  compared  with  Mie  theory  ( — )  for  a  =  1 
and  m  =  1,723  -  i  0.007. 
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Figure  a34.  Experimental  phase  function  data  (*)  for  a  particle  from 
138  for  ^  =  457.9  nm. 


A35.  Experimental  phase  function  data  (*)  for  a  particle  from  Run 


Figure  A37.  Experimental  phase  function  data  (*)  for  a  particle  from 
138  for  ^  =  457.9  nm. 
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Experimental  phase  function  data  (*)  for  a  particle  from  Run 
138  for  ^  =  488.0  nm. 


Figure  a39.  Experimental  phase  function  data  {*)  for  a  particle  from  Run 
138  for  A  =  514.5  nm. 


Figure  A40.  Experimental  phase  function  data  (*)  for  a  particle  from  Run 
147  for  X  =  457.9  nm. 


ANCLE  IN  DEGREES 

Experimental  phase  function  data  (*)  for  a  particle  from  Run 


Experimental  phase  function  data  (*)  for  a  particle  from  Run 
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